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The largest sample ever recorded of charged-current quasi-elastic (CCQE, v^+p — > u + +n) can- 
didate events is used to produce the minimally model-dependent, flux-integrated double-differential 

,2 _ 

cross section dT dcose ^ or v ^ incident on mineral oil. This measurement exploits the unprecedented 
statistics of the MiniBooNE anti-neutrino mode sample and provides the most complete information 
of this process to date. Also given to facilitate historical comparisons are the flux-unfolded total 
cross section a (E u ) and single-differential cross section on both mineral oil and on carbon by 
subtracting the Hp CCQE events on hydrogen. The observed cross section is somewhat higher than 
the predicted cross section from a model assuming independently-acting nucleons in carbon with 
canonical form factor values. The shape of the data are also discrepant with this model. These 
results have implications for intra-nuclear processes and can help constrain signal and background 
processes for future neutrino oscillation measurements. 



I. INTRODUCTION 

With the recent determination of the non-zero value 
of $i3 [1-4], present and future neutrino oscillation ex- 
periments will focus on measurements of the neutrino 
mass ordering and searches for leptonic CP violation. To 
reach discovery-level sensitivity to each of these effects, 
GeV-rangc v e and v e appearance must be observed in a 
long-baseline program with few-percent precision [5-9]. 
To facilitate such an ambitious program, the cross sec- 
tion for signal and background , v e , and v e charged- 
current processes must be known to high precision. This 
goal is commonly met by using a near detector to di- 
rectly measure the rate and shape of the unoscillatcd 
spectrum. However, if the cross sections are not inde- 
pendently and precisely understood, the extracted infor- 



mation may be model dependent and significantly biased 
from their true value [10, 11]. In the absence of a near de- 
tector, detailed knowledge of the contributing reactions 
is even more critical to the successful execution of these 
measurements. While the experimental and theoretical 
knowledge of GeV-range neutrino interactions on nuclear 
targets is improving, the experimental precision of inter- 
actions in this range is not better than 10% [12]. Of even 
more concern, as will be discussed in more detail, the fun- 
damental processes contributing to neutrino interactions 
with nuclear matter are not well understood. 

In general, anti-neutrino cross sections in the few-GcV 
region are not as well known as their neutrino coun- 
terparts, and in particular there are no charged-current 
anti-neutrino cross-section measurements below 1 GeV. 
In this work we present the first measurement of the anti- 
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neutrino charged-current quasi-elastic (CCQE) double- 
differential cross section with respect to kinematic prop- 
erties of the outgoing muon. These data are obtained 
using a muon anti-neutrino beam with mean energy 
(E p ) = 665 MeV and an exposure of 10.1 x 10 20 pro- 
tons on target (POT). This measurement represents an 
important step towards reaching the level of knowledge 
required for next-generation oscillation measurements. 

Apart from the valuable constraint these results pro- 
vide for future experiments seeking to use anti-neutrino 
events to measure the few remaining unknown funda- 
mental properties of neutrinos, the interpretation of the 
data will offer insight into an emerging puzzle. These re- 
sults significantly contribute to the body of experimen- 
tal information that suggest that the canonically-used 
model in neutrino generators of the Rclativistic Fermi 
Gas [13] (RFG) is insufficient for describing neutrino in- 
teractions in nuclear media. This model assumes the 
impulse approximation, under which nucleons housed in 
dense material are treated as quasi-free, independently- 
acting participants subject to a global binding energy 
and Fermi motion, while the surrounding environment 
behaves as simple spectators. In this formalism the in- 
teraction is parametrized by a set of tensor, vector and 
axial- vector form factors [14]. The vector form factors 
are measured in electron scattering data [15] while the 
axial-vector form factor is left to be empirically deter- 
mined by neutrino experiments and is typically assumed 
to take a dipolc form: 



F A = 



9A 



(1) 



where qa is measured from nuclear beta decay [16], Q 2 
is the squared four-momentum transfer and, while con- 
straints exist from pion electroproduction data [17], neu- 
trino experiments usually treat the axial mass Ma as a 
free parameter. 

By measuring the total rate of CCQE interactions 
and fitting the inferred Q 2 distribution, a variety of 
experiments employing bubble chamber detectors hous- 
ing mostly light nuclear targets typically produced con- 
sistent measurements of Ma- From these data, the 
averaged value is M A = 1.026 ± 0.021 GeV [17, 18]. 
With the discovery of neutrino oscillations, the use of 
light nuclear targets for the detection medium became 
impractical, as the statistics required to make high- 
precision oscillation measurements are much more eas- 
ily obtained using dense targets. With these relatively 
heavy nuclei and higher-precision detectors, more re- 
cent experiments have extracted values of Ma systemati- 
cally higher than 1.026 GeV [19-22]. Adding complexity, 
the modern heavy nuclear target experiment NOMAD 
has measured a value of Ma consistent with the bub- 
ble chamber analyses [23], while preliminary shape re- 
sults from the MINERVA experiment seem to also favor 
M A ~ 1 GeV [24]. 



An essential first step to understanding this apparent 
discrepancy is to recognize the particulars of the model 
dependence introduced by comparing values of Ma be- 
tween the many experiments. Important experimental 
differences that may contribute to the discrepancy in- 
clude disparate neutrino spectra, different neutrino de- 
tection technologies and the size of the nuclear media 
employed. Among the liberties taken to compare Ma val- 
ues across these scattering experiments include the dipolc 
form of Fa, various expectations of hadronic activity con- 
sistent with single-nucleon ejection and the previously- 
mentioned independent nucleon assumption implicit in 
both the formalism and in the inference of the Q 2 distri- 
bution. A possible reconciliation between the data sets 
has been proposed by offering a mechanism resulting in 
intra-nuclcar correlations of greater strength than pre- 
viously expected (see Ref. [25] and references therein). 
Such a mechanism is consistent with electron scatter- 
ing data [26, 27]. If this process is confirmed for weak 
interactions via neutrino scattering, its detailed under- 
standing will significantly expand knowledge of intra- 
nuclear behavior, and some neutrino oscillation results 
may need to be revisited [10, 11]. The best chance to 
definitively resolve this crucial ambiguity lies in the com- 
munity's ability and willingness to produce and com- 
pare model-independent information in both the leptonic 
and hadronic interaction sectors between experimental 
data and theoretical calculations. For this reason, the 
main result of this work is the double-differential CCQE 

cross section ( d , - a ] on mineral oil, where no as- 

sumptions about the underlying process is necessary for 
it's reconstruction. Regardless of the fundamental in- 
teractions contributing to the sample studied, this work 
reports the first cross-section measurements of ~ GeV 
anti-neutrinos and thus significantly advances the com- 
munity's preparedness to search for CP violation with 
neutrinos. 

This paper is organized as follows: the MiniBooNE 
experiment is described in Section II while Section III 
describes the model for neutrino interactions. The anal- 
ysis is presented in Section IV, and the conclusions are 
summarized in Section V. Appendix A presents a mea- 
surement of the charged-current background to the 
analysis sample, which exploits [i~ nuclear capture. Var- 
ious model-dependent CCQE cross sections are pro- 
vided in Appendices B and C, and Appendix D tabulates 
all cross-section results. 



II. THE MINIBOONE EXPERIMENT 

A. Beamline and flux 

MiniBooNE observes an on-axis neutrino flux from the 
Fermilab Booster Neutrino Beamline (BNB). A beam 
of 8.9 GcV/c momentum protons is extracted from the 
Booster synchrotron in bunches of 5 x 10 12 protons over 
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1.6 fj,s at an average rate of up to 5 Hz. A lattice of 
alternatively focusing and defocusing quadrupole mag- 
nets steers the proton spills into a beryllium target 71 
cm (1.75 interaction lengths) long. The protons collide 
with the target to create a spray of secondary particles. 
An aluminum electromagnetic horn surrounding the tar- 
get is pulsed to coincide with the p-Bc collisions, creating 
a toroidal magnetic field to focus mesons of the desired 
charge. For the data used in this analysis, the polarity 
of the magnetic horn is set such that negatively charged 
secondary particles are focused while those with positive 
charge are defocused. The accepted mesons are allowed 
to decay in a 50 m long air-filled hall, which terminates at 
a steel beam dump. The dominant decay modes of these 
mesons, mostly pions, produce muon neutrinos and anti- 
neutrinos. 

At MiniBooNE's request, the HARP experiment mea- 
sured pion production cross sections with a 8.9 GeV/c 
momentum proton beam on a 5% interaction length 
replica MiniBooNE target [28]. The HARP double- 
differential cross section in pion energy and angle min- 
imizes the model dependence of the BNB neutrino flux 
calculation [29] . A GEANT4-based model [30] takes these 
data as input and is used to predict the flux of neutrino 
and anti-neutrinos observed by the detector. The simu- 
lation considers proton transport to the target, p-Be in- 
teractions in the target including secondary interactions, 
meson production and their propagation through the 
magnetic field inside the horn, meson decay, and finally 
neutrino propagation to the detector. The uncertainty of 
primary tt~ production at the target is based exclusively 
on the HARP tt~ double-differential cross section [28]. 
Though the beryllium target used to collect the HARP 
data is substantially shorter compared to the MiniBooNE 
target (5% vs. 170% interaction lengths, respectively), 
the difference in ir production arising from the thickness 
between the two targets is calculated to be small. For 
the proton energies used by the BNB, roughly 90% of 
the neutrino beam is expected to come from the decay of 
primary 7r's [31], making the MiniBooNE flux prediction 
minimally dependent on the model for re-interactions in 
the long target. The anti- neutrino mode beam intersect- 
ing the detector is composed of 83.7% v^, 15.7% v^, 0.4% 
v e , and 0.2% v e . The and P M flux predictions are pre- 
sented in Figure 1. The electron- type neutrinos are ir- 
relevant to this analysis, but as the MiniBooNE detector 
is unmagnetized, the contribution represents a sig- 
nificant background. Furthermore, Figure 2 shows that, 
while the majority of v^s produced by ir~ decay are con- 
strained by the HARP measurement, most of the orig- 
inating from 7r + decay arise from a region not reported 
by HARP. In the analysis, the accepted flux of in the 
anti-neutrino mode beam is thus constrained using the 
observed rate of z/ M events in the MiniBooNE detector, 
as presented in Ref. [32] and Appendix A. These anal- 
yses constrain the knowledge of the flux to roughly 
15%. The fractional uncertainty of the flux predic- 
tion is around 7% at the interaction peak, roughly due in 



roughly equal amounts to errors on 7r production and 
the model that connects their production to the Vn flux. 




FIG. 1: (Color online) The MiniBooNE and flux predic- 
tions for anti-neutrino mode for the 10.1 xlO 20 POT exposure 



used in this analysis, 
vided in Table XI. 



Numerical values for the flux is pro- 
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FIG. 2: (Color online) Predicted angular distributions of pi- 
ons with respect to the incident proton beam (6^) producing 
Vfi and i> M events in the MiniBooNE detector in the anti- 
neutrino mode beam configuration. The v^/v^ event frac- 
tion is significantly larger than the flux fraction due to the 
respective cross sections. Distributions are normalized to 
10.1 x 10 20 POT, and arrows indicate the region of HARP 
data [28] constraints. Figure taken from Ref. [32]. 

Further details of the beamline and flux prediction are 
given in Ref. [29]. 



B. Detector 

The detector is a 12.2 m diameter sphere filled with 
818 tons of undoped mineral oil. The tank is optically 
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segregated into an inner signal region of radius 575 m 
and an outer veto shell of 35 cm thickness. Light pro- 
duced in the detector is collected by 1520 8-inch Hama- 
matsu photomultiplier tubes (PMTs), 1280 of which face 
into the signal region (11.3% coverage) while 240 are in- 
side the outer shell. Low activity in the veto region is 
required in physics analyses to ensure containment of 
charged particles produced by beam-induced neutrinos 
while also eliminating contamination from charged par- 
ticles entering the tank. 

Kept at ~ 20 °C, the mineral oil has a density of 
0.845 g/cm 3 with an index of refraction of 1.47. Under 
these conditions, charged particles with velocity f3 > 0.68 
produce Cherenkov radiation. Lepton particle identifica- 
tion and reconstruction is principally obtained through 
the pattern and timing of this prompt Cherenkov light. 
The PMTs have a quantum efficiency of ~ 20% and a 
timing resolution of ~ 2 ns, and the prompt Cherenkov 
component is easily separable from the delayed scintil- 
lation light present due to impurities in the oil. Four 
dispersion flasks at various locations in the detector are 
used to illuminate the signal-region PMTs with light from 
a pulsed laser. The laser data provides a calibration of 
PMT responses and allows an in situ measurement of 
light scattering properties over time. Throughout over 10 
years of MiniBooNE running, the observed energy scale 
has been stable to within 1%. 

PMT charge and time information is collected for a 
total of 19.2 /j,s beginning ~ 5 fis before the 1.6 fj,s long 
proton spill from the BNB. Cosmic ray muons stopped in 
the signal region prior to the start of the DAQ window 
may decay in time with the BNB spill, so PMT activity 
5 /is before proton delivery is monitored and used to mini- 
mize this contamination. Activity is recorded subsequent 
to the beam window for more than 10 fis to observe elec- 
trons from the at-rest decay of muons (hereafter referred 
to as "Michel" electrons) produced directly or indirectly 
through the primary neutrino interaction. 

The detector response to muons is calibrated using a 
dedicated system that independently measures the en- 
ergy and direction of cosmic ray muons up to 800 MeV. 
A scintillator hodoscope directly above the detector and 
seven scintillator cubes at various depths within the de- 
tector are used to track these particles. Each cube is 
connected by optical fiber to a PMT for readout. Sig- 
nals generated in the hodoscope and PMTs consistent 
with a muon stopping in a scintillator cube afford a di- 
rect calibration of the detector response to the range 
of muon kinematics most important to this analysis. 
These signals are used to verify muon reconstruction al- 
gorithms. Full details on the reconstruction is available 
in Ref. [33], while the detector is described in further 
detail in Rcf. [34]. 



III. PREDICTED NEUTRINO INTERACTIONS 

MiniBooNE uses the nuance [35] event generator 
to simulate neutrino interactions, nuance includes a 
comprehensive neutrino and anti-neutrino cross section 
model which considers known interactions in the neu- 
trino and anti- neutrino energy range from ~ 100 MeV to 
1 TeV. Ninety-nine reactions are modeled separately and 
combined with nuclear models describing bound nucleon 
states and final-state interactions to predict event rates 
and kinematics. 

Bound nucleons in the detector medium are described 
by the RFG [13]. This assumes the nucleons to be inde- 
pendent and quasi-free. Also specified is a hard cut-off 
in available struck nucleon energies as dictated by the 
exclusion principle. 

The neutrino interaction types most relevant to the 
current analysis are charged-current quasi-elastic (Sec- 
tion III A) and single pion production (Section IIIB). 
The neutrino-induced absolute cross sections for both 
processes have been measured at MiniBooNE using a flux 
prediction well determined by HARP data. These cross- 
section measurements are utilized in the anti-neutrino 
mode simulation. More broadly, to minimize the model 
dependence of the extracted CCQE cross section, each 
clear opportunity to constrain the backgrounds using 
MiniBooNE data was exploited. 

A. Charged current quasi-elastic scattering 

CCQE interactions are the most prevalent channel in 
MiniBooNE's energy range, and are predicted to account 
for ~ 40% of all events. Their simulation in this analysis 
is chosen based on results from the MiniBooNE and 
CCQE data sets. The formalism is described by a 
relativistic Fermi Gas model [13], and with a few empiri- 
cal parameter adjustments, this model adequately repro- 
duces the kinematics of both CCQE data sets [19, 36]. 
Through the procedure to correct for detector resolution 
effects, this choice only mildly affects the shape of the 
extracted true muon kinematics, while the normalization 
of the distribution is entirely unaffected. It will be shown 
later that this effect is negligible compared to other sys- 
tematic uncertainties. 

The vector and tensor components of the interaction 
are constrained by data from electron scattering experi- 
ments and a non-dipole form is taken based on the results 
of Rcf. [15]. As shown in Equation 1, the axial- vector 
form factor assumes a dipole form and contains the em- 
pirical "axial mass" parameter Ma- In this analysis, the 
value of Ma is chosen based on results from neutrino 
interactions. 

As CCQE interactions exclusively interact with nu- 
cleons bound in carbon, M A S ' C = 1.35 GeV together with 
a Pauli blocking adjustment, k = 1.007 is sufficient to de- 
scribe the kinematics of such events based on a fit to the 
MiniBooNE data [19]. With the kinematics character- 
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ized by this adjusted prediction, the total cross section 
for the simulated process is subsequently corrected to the 
observed normalization in data. In this way, the details 
of the observed CCQE data are reproduced in the 
present simulation for this process. 

The MiniBooNE mineral oil is composed of C n H 2n +2 5 
n ~ 20, and so CCQE scattering off of both bound and 
quasi- free protons are accessible to £v's. For the hydro- 
gen scattering component, M^ ff,H = 1.02 GeV is cho- 
sen based on the body of experimental results for the 
CCQE process incident on light nuclear targets [17, 18]. 
For interactions with protons bound in a carbon nucleus, 
the binding energy (Fermi momentum) is set to 30 MeV 
(220 MeV) based on electron scattering data for the QE 
process [37]. Along with the same CCQE model param- 
eters measured in the data of AI A B ' C = 1.35 GeV and 
K = 1.007, these choices are adopted for CCQE in- 
teractions on carbon. This choice is made exclusively 
due to observed agreement between this model and the 
MiniBooNE CCQE data [36]. Note that, due to the 
axial- vector interference term in the formalism, the kine- 
matics of CCQE feature a softer momentum-transfer 
spectrum and so, in the RFG, the same value of k has a 
larger effect on CCQE compared to CCQE. More 
importantly, it will be shown later that the extracted £y, 
CCQE double-differential cross section is only negligibly 
affected by these choices. 

The superscript "eff", short for "effective", on Ma is 
used throughout this work to allow for the possibility that 
nuclear effects are responsible for the apparent discrep- 
ancy between the light target results and those from the 
more recent experiments using dense nuclear material. 
As discussed in the introduction, this is also motivated 
by theoretical work that predicts an extra class of events 
whose contribution to the CCQE sample in Cherenkov 
detectors, such as MiniBooNE, enters due to the lack 
of requirement on hadronic activity [25]. The letter fol- 
lowing the "eff' identifies the relevant nucleon target in 
MiniBooNE's hydrocarbon medium: H for the quasi-free 
hydrogen targets, and C for those bound in carbon. 

B. Pion production 

The majority of single pion production (yi + TV — > 
I + 7T + TV') events at MiniBooNE energies are medi- 
ated by baryonic resonances. The formalism to describe 
these events is taken from the Rcin-Schgal model [38], 
where the rclativistic harmonic oscillator quark model is 
assumed [39]. The production of A(1232) is dominant in 
the energy range spanned by MiniBooNE, but 17 other 
and higher-mass resonances are also considered. 

The charged-current single pion channels for (y^ + 
TV ->■ n~ + 7T+ + TV, "CC17T+") and {v^ + TV -> 
fi + + tt~ + TV, "CC17T - ") dominate the pion-producing 
interactions contributing to the CCQE sample. The 
CCl7r + events enter from the content of the beam. 
The CCl7r~ background results from interactions, 



and their presence in the CCQE sample is mostly due 
to stopped 7T _ capture in the nuclear medium. Stopped 
ir~ capture in the presence of carbon is ~ 100%, so they 
are not separable from the CCQE sample through ob- 
servation of 7T~ decay. In the current analysis, the Rein- 
Schgal prediction for both classes of events is adjusted to 
reproduce the kinematic distributions measured in Mini- 
BooNE neutrino-mode CC1tt + data [19, 40]. 



C. Final state interactions 

Subsequent to a neutrino interaction involving a nu- 
cleon bound in carbon, nuance propagates the outgo- 
ing hadrons including nuclcons, mesons and baryonic 
resonances, and simulates their re-interaction as they 
exit the nucleus. The initial interaction model employs 
the impulse approximation which assumes an instan- 
taneous exchange with independent nucleons. Subse- 
quent to the initial neutrino or anti-neutrino interac- 
tion, particles produced inside the nucleus are propa- 
gated step-wise in 0.3 fm increments until they emerge 
from the ~ 2.5 fm radius sphere. Intermittently, the 
probability for hadronic rc-intcraction is calculated us- 
ing a radially-dependent nucleon density distribution [41] 
along with external 7r — TV, TV — TV cross-section measure- 
ments [42]. For A re-interactions (A + TV ->• TV + TV), an 
energy-independent probability of 20% (10%) is taken 
for A+ + TV, A + TV (A++ + TV, A" + TV) based on K2K 
data [43] and is assigned 100% uncertainty. The dom- 
inant final-state interactions affecting this analysis are 
pion charge exchange (7r ± + X <H> n° + X ) and absorp- 
tion (tt ± + X -> X). 



IV. ANALYSIS 

This section describes the extraction of the CCQE 
double-differential cross section. It is necessary to first 
identify the experimental complications that distinguish 
this measurement from the MiniBooNE CCQE result. 

Though the same detector, reconstruction and event 
selection are used for the [19] and CCQE analy- 
ses, subtleties related to the detector material and the 
different beam configuration result in substantially dif- 
ferent sample content in both the signal and background 
processes. Due to leading-particle effects at the beryl- 
lium target, the mean energy of the flux in anti- 
neutrino mode (shown in Figure 1) is appreciably lower 
((Ep) = 665 MeV) compared to the flux in neutrino 
mode {{E v ) = 788 MeV). The content of the two CCQE 
signal samples are also fundamentally different since 
CCQE events arise from interactions with protons while 
CCQE events involve interactions on neutrons. The 
hydrocarbon nature of the detection medium provides a 
mix of bound and quasi-free interaction targets for 
CCQE, while CCQE involves only bound nuclcons. 
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The two interaction types for CCQE are not sepa- 
rable, and so the sum of all za, CCQE interactions are 
treated as the signal for this analysis. However, as his- 
torical data on mostly light targets are adequately de- 
scribed with Ma ~ 1 GeV, results evoking this model 
to subtract the quasi-free CCQE content are given in 
Appendix B. 

Backgrounds in this analysis also offer unique compli- 
cations, as mentioned in Section II A and expanded in 
the next section. Broadly, the analysis sample is formed 
with a simple selection that requires the prompt muon 
be contained in the detector and that its decay is ob- 
served. The dominant backgrounds with this selection 
are charged-current (CC) and CCl7r~ interactions. 
The ffj, CCQE contribution is indistinguishable event-by- 
event from v^, CCQE; however statistical measurements 
of their overall rate and shape discussed in Section IV A 
constrain the knowledge of this background to ~ 15%. 
The CCl7r~ contamination enters from the capture pro- 
cess on carbon nuclei and is known less precisely, as it is 
not separable in the data from CCQE. Furthermore, 
there are no measurements of this process in external 
data sets at the MiniBoNE energy range. Due to the size 
of the CCQE and single-pion backgrounds, the signal 
purity is only 61% in this work, compared to 77% for 
the i/fj, CCQE analysis. Multiple dedicated analyses and 
comparisons were necessary to reduce the uncertainty on 
these processes to a manageable level, and result, 
the final uncertainty on the extracted cross sections 
are dominated by the level of P M flux uncertainty. 



A. Constraints on background processes 



rated. These techniques include the use of [i~ nuclear 
capture, 7r~ nuclear capture, and angular differences be- 
tween and in CCQE interactions. The analysis 
based on \x~ capture is described in Appendix A and 
the other analyses are presented in Ref. [32]. The \i~ 
capture analysis takes advantage of the fact that ~ 8% 
of ^-induced CC interactions on carbon that do not 
lead to a decay electron. Nuclear capture of tt~ also 
affords sensitivity to the beam content. The second- 
most prevalent interaction type in the MiniBooNE de- 
tector is CC single-pion production, which produces a 
7r + in the case of scattering and a 7r~ for reac- 
tions. As almost all stopped ir~ are absorbed in the hy- 
drocarbon medium [44] , the sample consisting of a single 
muon and two decay electrons (one each from the prompt 
muon and the pion decay chain) is predominantly due to 
events. Finally, the observed angular distribution of 
CCQE events is fit to a combination of and £y, events, 
where interactions are predicted to be much more 
forward-going with respect to the beam direction. This 
last analysis is dependent on details of the RFG predic- 
tion for CCQE scattering, and its results are not used 
in the subtraction of the background. In the future, 
where CCQE and CCQE-like interactions should be bet- 
ter understood, this technique could provide a valuable 
constraint. The results from these analyses are summa- 
rized in Figure 3, where the nominal and highly-uncertain 
prediction of the flux in the anti-neutrino mode beam 
appears to be roughly 20% high in normalization, while 
the energy dependence seems to be well modeled. Based 
on the results from the fi~ and tt~ nuclear capture anal- 
yses, the Vy, flux in the anti-ncutrino mode beam is cor- 
rected by a scale of 0.77 with an uncertainty of 0.13. 



The largest background in the CCQE sample is the 
Vy contamination. Moreover, as Figure 2 shows, the ma- 
jority of 7r + particles contributing to the beam are pro- 
duced at small angles with respect to the incoming pro- 
tons (and so are affected less by the magnetic field) and 
thus their contribution to the anti-neutrino mode beam 
is mostly unconstrained by the HARP hadroproduction 
data. As MiniBooNE is non-magnetized, this motivated 
a dedicated study of the beam content using statis- 
tical methods. Three techniques, described in detail in 
Ref. [32] and Appendix A, were used to measure this cru- 
cial background for the MiniBooNE data. Two of these 
measurements are largely model independent, and the fi- 
nal fractional uncertainty on the beam contribution 
is ~ 15%. These analyses are the first of their kind and 
their uncertainty has reduced the error on the CCQE 
cross section due to z/ M interactions to a sub-dominant 
uncertainty. 

The three measurements of the flux of the anti- 
neutrino mode beam exploit various differences between 
charged-current and v^. processes to statistically mea- 
sure their respective contributions. Broadly, these mea- 
surements are executed by performing rate analyses on 
samples with the and content statistically sepa- 



1.2 — 



1 — 




- • CC1jc + 
0.2 ■ CCQE 
Z ▼ H" capture 

, I . I . . I . . , . . .1 , I . 
0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 ALL 
generated E v (GeV) 

FIG. 3: (Color online) Summary of the results from three 
techniques used to measure the v^. flux in the anti-neutrino 
mode beam. The "/x _ capture" analysis is described in Ap- 
pendix A, and the others in Ref. [32]. The measurements are 
given relative to an extrapolation of HARP data into a region 
where no relevant hadroproduction data exists (shown as a 
dotted line at unity). 
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Chargcd-current single tt~ production constitute the 
second-largest background to this analysis, accounting 
for ~ 15% of the sample. These events enter through 
a different and experimentally disadvantageous mech- 
anism compared to the analogous process for the 
CCQE sample. Single-pion events induced by typi- 
cally give rise to Michel electrons through the decay chain 
7r + — > /i + — > e + of stopped pions, which can be observed 
and used to reject these events. However, an appreciable 
number of 7r + are destroyed in flight through the nuclear 
absorption process (ir + + X — s- X'). In contrast, almost 
all single-pion events from interactions enter the 
CCQE sample since ir~ that are not absorbed in flight 
stop and undergo nuclear capture with ~100% efficiency 
While this fortuitously allows for the CCl7r + -based 
flux measurement, this also implies CC17T - events arc 
not separable from the CCQE sample. This is in con- 
trast to the MiniBooNE CCQE analysis, where the 
single-pion events tagged through the observation of an 
additional Michel allows a direct constraint of the rate 
and kinematics of CCl7r + events. A correction was thus 
measured and applied to the background prediction for 
the Vfj, CCQE sample [19]. This constraint is applied 
by correcting the CCl7r + events according to the ob- 
served reconstructed momentum transfer. The correc- 
tion is shown in Fig. 7b of Ref. [19]. In the absence 
of such a measurement for CCl7r _ interactions, the con- 
straint obtained in neutrino mode for CCl7r + is ap- 
plied to the CCl7r~ Rein-Sehgal prediction described in 
Section IIIB. Figure 4 shows this prediction agrees well 
with an external calculation [45] for such events. This al- 
ternate model is implemented in NUANCE and is based on 
an extension of the Rein-Sehgal model found in Ref. [46] . 
This updated calculation includes muon mass terms and 
a modified vector form factor to yield better agreement 
with world pion production data. Consistency between 
these two predictions for CCl7r _ production, along with 
the level of agreement between the extended Rein-Sehgal 
calculation and the MiniBooNE CCl7r + data (shown in 
Ref. [45]) suggests an uncertainty of 20% is sufficient for 
the CC17T - background. 

Based on results from the neutrino-mode CCl7r° 
analysis [47], the small contribution from CCl7r° events 
induced by both and are increased by a factor of 
1.6 relative to the NUANCE prediction. The generated 
predictions for all other interactions, accounting for < 3% 
of the sample, are not adjusted. 



B. Reconstruction and analysis sample 

The identification of CCQE candidate events re- 
lies solely on the observation of a single muon and no 
final-state pions. Muon kinematics are obtained by the 
pattern, timing, and total charge of prompt Cherenkov 
radiation collected by PMTs. A likelihood function oper- 
ating under a muon hypothesis is compared to the topol- 
ogy and timing of the observed PMT hits. This likcli- 



CC17T events in 7 CCQE sample 




FIG. 4: (Color online) Three calculations for the CCl7r~ 
background contribution to the CCQE sample as a func- 
tion of cos# M , where # M is the angle of the muon relative 
to the neutrino direction. The "central value" distribution 
corresponds to the nominal Rein-Sehgal [38] prediction for 
CCl7r~ events in MiniBooNE constrained by the observed 
kinematics in the neutrino-mode CCl7r + sample. This 
agrees well with a more recent calculation ("Berger-Sehgal 
extended model" ) that is based on an improved version of the 
Rein-Sehgal model. For comparison, the bare Rein-Sehgal 
prediction for CCl7r~ events is also shown. Distributions are 
normalized to 10.1 x 10 20 POT. 



hood function predicts hit patterns and timing based on 
the interaction vertex and the momentum four-vector of 
the muon. The likelihood function simultaneously varies 
these seven parameters while comparing to the observed 
PMT hits. The parameters from the maximized likeli- 
hood function yield the reconstructed muon kinematics. 
Integrated over the spectrum of observed muons, the res- 
olution of this reconstruction for muon energy (angle) 
is roughly 8% (2 degrees) [33]. The direct and high- 
resolution observation of muon properties motivates the 
choice to present the CCQE cross section as a func- 
tion of muon kinematics as the main result of this work, 
while the unprecedented statistics of the data set ana- 
lyzed yields sensitivity to previously unprobed regions of 
the interaction. 

As in the CCQE work, no requirement is made on 
hadronic activity. This is an important distinction from 
the CCQE definitions used by other experiments [21, 23], 
where a single proton track may be required for CCQE 
selection. However, note that in the case of CCQE 
scattering, where a single ejected neutron is expected, 
the experimental definition used by tracking detectors 
is largely based on a single muon track. Therefore, in 
general, the selection used by tracking detectors and 
Cherenkov-based measurements such as this one for 
CCQE follow each other more closely as compared to the 
Vfj, case. 

Event selection for this work is thus identical to that 
used in the MiniBooNE CCQE analysis [19]. Table I 
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TABLE I: Sample purity and detection efficiency for all CCQE events, which are due to a mix of scattering on bound and 
quasi-free nuclear targets. Efficiencies are normalized to events with a generated radius r < 550 cm. 



Cut # 


Description 


Purity 

\ /0 ) 


Efficiency 

\ /0 ) 





No cuts 


32.3 


100 


1 


Veto hits < 6, all subevents 


27.6 


50.8 


2 


First subevent: in beam window 4000 < T(ns) < 7000 


27.7 


50.3 


3 


First subevent: muon kinetic energy T M > 200 MeV 


36.9 


44.0 


4 


Only two subevents 


48.4 


38.8 


5 


First subevent: reconstructed vertex radius R < 500 cm 


49.2 


32.6 




Distance between subevent reconstructed vertices > 500 cm/GeV x T M - 100 cm 






6 


Distance between subevent reconstructed vertices > 100 cm 


54.3 


30.6 


7 


First subevent: log-likelihood (fj, / e) > 


61.0 


29.5 



provides cumulative purity and efficiency values for the 
selected sample. Sample selection is based on require- 
ments of temporally-correlated collections of PMT activ- 
ity (or PMT "hits") referred to as "subevents". A hit 
is any PMT pulse passing the discriminator threshold of 
~ 0.1 photoelectrons, and a cluster of at least 10 hits 
within a 200 ns window with individual hit times less 
than 10 ns apart defines a subevent. Two or fewer spac- 
ings between 10 - 20 ns among individual hit times are 
also allowed. The primary requirement to identify 
CCQE events is two and only two subevents, due domi- 
nantly to Cherenkov light from the prompt muon and its 
decay positron: 



/t + n 



(2) 



The difference in average PMT hit time between the 
two subevents is given in Figure 5 and shows both the 
characteristic lifetime of muons in the sample and the 
effect of the subevent definition on CCQE detection for 
quickly-decaying muons. The selection criteria are enu- 
merated in Table I. Cut 1 enforces containment of 
charged particles produced inside the detector while also 
rejecting incoming charged particles. Cut 2 requires the 
muon subevent be correlated with the BNB proton spill. 
Cut 3 ensures the first subevent is not a Michel electron 
and avoids a region of muon energy with relatively poor 
reconstruction. Cut 4 eliminates most neutral current 
events and rejects most interactions with final-state 7r + . 
Cut 5 further enhances the reliability of the reconstruc- 
tion by reducing sensitivity to PMT coverage. Cut 6 
ensures the measurements of the muon energy and the 
subevent vertices are consistent with the production and 
subsequent decay of a minimum ionizing particle. This 
cut rejects many events where the Michel is not associ- 
ated with the primary muon, mainly CCl7r + and NCl7r + 
events where the second subevent is a decay positron 
from the ir + decay chain. Cut 7 requires that the can- 
didate primary muon is better fit as a muon than as an 
electron. This cut reduces the background from most 
processes, most notably from CCl7r + and CCl7r~. 




1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 
Michel avg. time - ii avg. time (ns) 



FIG. 5: (Color online) Time separation in the average PMT 
hit time between the two subevents in the CCQE sample. The 
deviation from an exponential form at low time differences 
is due to the subevent separation requirement. Simulation 
is normalized to data, and statistical errors are shown with 
data. 



For this selection, the neutrino interaction assumptions 
detailed in Section III, the constrained backgrounds de- 
scribed in Section IV A and 10.1 x 10 20 POT, the sample 
consists of 71,176 events with CCQE purity (detection 
efficiency) of 60.3% (29.5%). Table II presents a sum- 
mary of the Vp, CCQE sample, while Figure 6 shows how 
these channels contribute to the reconstructed kinemat- 
ical distributions of the final state muon. 



C. Cross-section calculation and uncertainties 

The flux-integrated double-differential cross section 
per nucleon in the i th bin is given by: 



da 



J2 3 Uij (dj-bj) 



dT^ d [ cos M ] y . (AT,), (A [ cos a $ N 



, (3) 
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TABLE II: Summary of the CCQE sample. Contributions 
reflect all adjustments to simulation based on constraints from 
MiniBooNE data. 





j- „ „„„ j- _ J\ T) /'VI' 

integrated _r U 1 


1U.1 X 1U 




energy- integrated Dn flux 


z.yo X 11) / cm 




t'u w v.. V ^ I j LetllLllLLeiLfc! cvfcJIlLb 


71 1 7fi 




CCQE efficiency (R < 550 cm) 


29.5% 




Interaction channel 


Contribution (%) 




-(- p — ^ -(- ti (bound p) 




V 


^ + p — > + n (quasi-free p) 


1 ( .1 




Vfi -\- n + p 


16.6 


v» 


f AT — >• + JV + 7r~ (resonant) 


10.4 


Vy 


f AT — > |U~ + AT + 7r + (resonant) 


3.8 




+ A ^ p> + + A + n~ (coherent) 


3.3 




v„. + N ->■ n+ + N + tt° 


2.8 




P M + p -> a* + + A 






^ M + n — ► m + + 


2.0 




^ M + p -»■ p, + + E° 






All others 


0.8 



where dj refers to data, bj the background, Uij is an 
unfolding matrix connecting the reconstructed variable 
index j to the true index i, is the detection efficiency, 
AT M and A [cos M ] the respective bin widths, $ the inte- 
grated Vy exposure, and N the number of proton targets 
in the volume studied. 

The unfolding matrix Uij is based on the Bayesian 
method proposed in Ref. [48] to account for reconstruc- 
tion biases. The high-sensitivity resolution of the recon- 
struction used to identify muon kinematics leads to only 
mild corrections. However, this procedure does introduce 
some dependence on the generated muon kinematics of 
Vp CCQE interactions. This bias is evaluated by unfold- 
ing the data with 100 different versions of EZy generated 
using a conservative range of CCQE model parameters. 
The bias introduced by the Bayesian unfolding method 
for the cross sections reported here is found to be negli- 
gible. 

It is important to note that by directly subtracting the 
background from data in the reconstructed distribution, 
this cross-section extraction procedure is entirely inde- 
pendent of the normalization of the generated signal 
CCQE processes. That is, though the RFC with a large 
value for the effective axial mass is assumed by simula- 
tion, the extracted cross section is not affected by this 
choice. 

Systematic uncertainties are evaluated by forming an 
error matrix that propagates correlated uncertainties on 
parameters and processes that affect CCQE inter- 
actions onto the calculated cross section. The covari- 
ance matrix is constructed by first forming a distribu- 
tion of weights corresponding to simulated excursions set 
by Gaussian variations of parameters and measurements 
within their associated error. These weights are then 
used to recalculate the double-differential cross section in 



6000 



5000 



4000 



> 3000 



2000 



1000 




CCQE (p from 12 C) 
v, CCQE (p from HJ 
| | v M CCQE 
\ | CCtt 

|] CCti + 

j | all others 

(a) 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

MGeV) 




FIG. 6: (Color online) Predicted sample composition as a 
function of reconstructed muon kinematics in the CCQE 
sample. The top figure (a) shows the kinetic energy of the 
muon, while the bottom (b) shows cos 9^, where # M is the 
muon direction relative to the incoming neutrinos. Distribu- 
tions are normalized to the exposure of 10.1 x 10 20 POT. 



Eqn. 3, replacing the central-value MC valuations with 
the excursion values for terms appropriate to each sys- 
tematic uncertainty. The difference of these alternate 
cross-section calculations compared to the "best guess" 
distribution forms the covariance matrix: 



s=l 



Nf v ) x (m - Nf v ). (4) 



Here K simulation excursions are used, A^ s is the re- 
weighted cross section value corresponding to the s th sim- 
ulation set and N cv represents the simulation central 
value. This technique is further described in Ref. [49]. 
For uncertainties on processes with correlated errors, typ- 
ically K = 100 while K = 1 is sufficient for uncorrelated 
errors. Systematic uncertainties requiring correlated er- 
rors are the production of ir~ in the proton beam target, 
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the connection between n~ production and the focused 
9^ beam, optical transport in the detector, final-state in- 
teractions, and the bias due to the unfolding procedure. 

As mentioned in Section II A, the uncertainty on the 
production of parent 7r~ at the beryllium target is 
driven by the HARP data [28] and the absolute i> M flux 
prediction is minimally dependent on the hadroproduc- 
tion model. Subsequent to 7r~ production, errors on 
the processes that culminate in the beam include the 
amount of delivered POT, optics of the primary beam, 
magnetic focusing, and hadronic interactions in the tar- 
get and the enclosing horn. More details on uncertainties 
of the flux prediction are available in Ref. [29]. Uncer- 
tainties on the model for optical transport in the detec- 
tor are based on both external and in situ measurements 
of light attenuation, scintillation strength, and the re- 
fractive index of the oil [50]. For this uncertainty, 70 
samples generated with variations of 35 parameters that 
describe the optical model are used to find the uncer- 
tainty propagated to the measurement. The most im- 
portant final-state interactions affecting the composition 
of the Vp CCQE sample are the pion charge-exchange 
(7r ± + X o 7T° + X') and absorption (^ + X ->■ X') 
processes. The uncertainty on pion charge- exchange (ab- 
sorption) inside the nucleus is set to 30% (25%) based on 
the difference between the nuance prediction and exter- 
nal data [51]. The intermedium processes are evaluated 
separately with 50% (35%) fractional uncertainty based 
on comparisons with the GCALOR prediction and the 
same external data. The final correlated systematic error 
evaluates the bias introduced by the Bayesian unfolding 
procedure, where 100 different matrices Uij are generated 
within M A = 1.35 ± 0.35 GcV and k = 1.007 ± 0.007. 
The negligible bias found when the data is extracted with 
these alternate matrices assuming a conservative range of 
CCQE parameter values assures this cross-section mea- 
surement is largely independent of the CCQE interaction 
model. 

Uncertainties described by a single excursion from 
the simulation central value include errors due to de- 
tector PMT response and on background processes not 
due to final-state interactions. Large sets of simulation 
are generated separately to evaluate rate biases due to 
uncertainties on phototube discriminator threshold and 
the correlation between pulse time and delivered charge. 
The background processes are grouped into three classes: 
CC17T - events, those induced by u^, and all non-CCQE, 
non-CC single-pion events. Note these groups are not 
mutually exclusive and all constraints are described in 
Section IV A. Based on consistency of the prediction us- 
ing an extrapolated CC17T constraint with a robust exter- 
nal model for the CCl7r~ background, these events are 
assigned 20% uncertainty. All background events are 
subject to the measured 15% fractional flux uncertainty. 
The cross sections for the CCQE and CCl7r + processes 
are directly measured by MiniBooNE data [19, 40], and 
so only their flux is uncertain. The small contribution 
from non-CCQE, non-CClTr 1 * 1 processes are largely un- 



constrained and assigned 30% cross-section uncertainty. 

The overall size of these covariance matrices can be ex- 
pressed with a single number, representing the normal- 
ization uncertainty of each error. Using the sum rule for 
variances and covariances, the total normalization uncer- 
tainty can be thought of as the error on the cross section 
if the measurement consisted of a single bin: 



5D T /D T 



\ 



(5) 



where Dt = represents the double-differential 

i 

cross section measurement summed over each kinematic 
region i. This is also commonly referred to as the uncer- 
tainty on the scale of the measurement. Table III shows 
the contributions of various errors to the total normal- 
ization uncertainty. 

The covariance matrix can also be used to separate 
the correlated normalization uncertainties from the to- 
tal error, leaving information related to how much the 
shape of the observed data may vary within the system- 
atic errors [52]. These uncertainties are identified by 
first defining a data vector V with entries correspond- 
ing to the observed relative normalization of each bin: 
Vi = {D 1 /D T ,D 2 /D T ,--- ,D n /D T ,D T }. Notice this 
vector has dimension n + 1, where n is the number of 
bins measured. The covariance matrix Q for this new 
vector V involves the Jacobian matrix of partial deriva- 
tives J and is given by: 



Qki = E J * M U ■>'.> = E Qfj-Mii JjJT ■ ( fi ) 



The diagonals of the matrix Q are related to the 
shape uncertainty in each kinematic bin. For entries 
{1,2,--- ,n}, 



N 2 

N 2 ^ 13 

1 ij 



(7) 



(<5-Dfc, shaped 



As the full covariance matrix M for the double- 
differential cross section is in principle a four-dimensional 
object with over 100,000 entries, the combination of the 
total normalization error and the bin-by-bin shape error 
is the preferable method to report the complete experi- 
mental uncertainty. This is argued more completely in 
Ref. [48], and Rcf. [53] provides an example of how do 
use this information in the context of a fit to these data. 

The main results of this work is the CCQE double- 
differential cross section on mineral oil. However, as the 
majority of the bubble-chamber CCQE analyses using 
light targets for the interaction medium are adequately 
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TABLE III: Normalization uncertainty for various sources of 
error for the CCQE cross section on mineral oil. 



Uncertainty type 


Normalization uncertainty (%) 


p flux 


9.6 


Detector 


3.9 


Unfolding 


0.5 


Statistics 


0.8 


Vfi background 


3.6 


CCl7r~ background 


4.0 


All backgrounds 


6.4 


Total 


12.9 



described with Ma ~ 1 GeV [17, 18], the cross section 
on carbon only is found by assuming this value to sub- 
tract the quasi-free hydrogen content of the CCQE 
data. This alternate cross section is calculated by includ- 
ing hydrogen CCQE events in the background term 
bj in Equation 3, while the other terms in the calculation 
based on the signal definition now are based on only 
CCQE events involving protons bound in carbon. Most 
notably, this reduces the number of interaction targets in 
the fiducial volume. 

Informed by the results of fits to the light-target CCQE 
experiments, Mf' U = 1.026 ± 0.021 GeV [17, 18] is as- 
sumed and subtracted from the data. Systematic error 
due to this background is evaluated with the method 
described earlier in this section with K = 100 throws 
against the 0.021 GeV uncertainty. Including this ad- 
ditional error and, more importantly, considering the 
lower sample purity for this alternate definition of signal 
events, the fractional normalization uncertainty increases 
to 17.2%. 



D. Results 

The CCQE double-differential flux-integrated cross 
section on mineral oil is shown with shape uncertainty 
in Figure 7 and it's one-dimensional projections are com- 
pared to RFC predictions in Figure 8. The configura- 
tion with the hydrogen content subtracted is given in 
Appendix D and may be more readily compared to the- 
oretical calculations for CCQE interaction on carbon, 
such as in Refs. [54-58]. Bins in the kinematic region 
-1 < cos 6^. < +1 and 0.2 < T M (GeV) < 2.0 
are reported if they meet the statistical requirement of 
at least 25 events in the reconstructed and background- 
subtracted data term (dj — bj) in Equation 3. If this 
threshold is not met, no measurement is reported. As no 
explicit assumptions about the underlying interaction are 
necessary to reconstruct muon kinematics, this result is 
nearly model independent. Since some background pro- 
cesses are not directly constrained by data, most notably 
CCl7r~, Appendix D tabulates the subtracted data. 




FIG. 7: (Color online) The P M CCQE double-differential cross 
section with shape uncertainty. Numerical values for this 
cross section and its uncertainty is provided in Tables XII 
and XIII, respectively. 



V. CONCLUSION 

This work presents the first measurement of the 
CCQE double-differential cross section in terms of muon 
angle and energy. This measurement is also the first 
charged-current cross-section measurement with the ma- 
jority of interactions with < 1 GeV. This cross section 
is the least mo del- dependent measurement possible with 
the MiniBooNE detector and is thus the main result of 
this work. 

It is clear in Figure 8 that the RFG model with 
canonical assumptions does not adequately describe these 
data neither in shape nor in normalization. Consis- 
tent with other recent CCQE measurements on nuclear 
material [19-22], a significant enhancement in the nor- 
malization that grows with decreasing muon scatter- 
ing angle is observed compared to the expectation with 
M A = 1.0 GeV. 

These data find tension with the NOMAD P M CCQE 
results, which are described both in shape and normaliza- 
tion by Ma = 1-06 ± 0.12 [23]. This tension is common 
among the CCQE analyses from the two experiments. 
However, care should be taken in comparing model- 
dependent results among experiments with such different 
neutrino fluxes and detector technologies. A definitive 
unification of these apparently discrepant data sets will 
require the continued increase of both experimental and 
theoretical activity surrounding this topic. Fortunately, 
many experiments at a variety of neutrino energies ca- 
pable of making high-resolution, mo del- independent neu- 
trino and anti-neutrino CCQE measurements with differ- 
ent detector technologies and nuclear media using both 
neutrino and anti-neutrino beams currently have data or 
will soon. These include MINERVA [59], SciBooNE [60], 
MicroBooNE [61], ArgoNeuT [62], ICARUS [63] and the 
T2K [4] and NO^A [5] near detectors. 
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Finally, a novel and crucial evaluation of the back- 
ground in this work is presented in Appendix A. In 
the absence of a magnetic field, this analysis and those 
described in Rcf. [32] measure the flux of the anti- 
neutrino mode beam with ~ 15% fractional uncertainty. 
These techniques could be used in current and future 
neutrino oscillation programs, particularly when mod- 
est charge identification is sufficient to meet the physics 



goals [64]. 
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FIG. 8: (Color online) Projections of the double-differential 
cross section in muon kinetic energy T M for various scattering 
angles cos(9 M . As indicated, both RFG predictions assume an 
effective axial mass of 1.02 GeV for the quasi-free hydrogen 
component of the data, while two choices for CCQE model 
parameters relevant to nuclear scattering are shown. Data 
with cos 9^ < -0.6 have insufficient statistics and this region 
is not shown. Numerical values for this cross section and its 
uncertainty is provided in Tables XII and XIII, respectively. 
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Appendix A: Measurement of flux in 
anti-neutrino mode using pT capture 

1. Introduction 

MiniBooNE uses dedicated hadroproduction measure- 
ments from the HARP experiment [28] to predict the 
and fluxes for the anti-neutrino mode beam. However, 
as shown in Figure 2 in Section II A, most of the flux 
arises from the very forward-going region of ir + produc- 
tion and is not well-constrained by the HARP measure- 
ments. Fortunately there are several ways to determine 
the content of the beam directly from MiniBooNE 
data. Two such analyses are described in Ref. [32] and a 
third is presented in this appendix. These analyses show 
that, in the absence of a magnetic field, the and 
content can still be modestly separated using statistical 
methods. 

The measurement of the flux in the anti-neutrino 
mode beam described in this appendix exploits the asym- 
metry in the production of decay electrons between yT 
and fi + in nuclear material. The results are consistent 
with and complementary to those of Ref. [32] . 



2. Muon capture model and event selection 

The model for fi~ capture and the processes that can 
obscure it's rate in the MiniBooNE detector is described 
in this section, followed by details on the analysis sam- 
ples studied. In mineral oil, stopped /i - arc captured on 
carbon nuclei with a probability of (7.78 ± 0.07)% [65]. 
In such capture events, typically little or no extra activ- 
ity is observed in the detector. However, the low-energy 
neutron and photons from the primary capture reaction 
as well as de-excitations of the boron isotope may be 
energetic enough to produce a Michel-like event. The 
simulated production of these particles is based on the 
measurements of Rcfs. [66-71], and the model that propa- 
gates these particles and possible re-interactions through 
the MiniBooNE detector estimates 6.60% of /i~ capture 
events lead to activity similar to a low-energy Michel. 
Thus, the apparent \i~ nuclear capture probability in the 
detector is predicted to be 7.78 x (1 - 6.60%) = 7.26 ± 
0.20%, where the uncertainty is substantially increased 
to recognize the model dependence of the rate to regain 
Michel-like events following /i~ capture. This rate is par- 
tially constrained by the calibration procedure described 
in Section A 3, and it will be shown that the assigned un- 
certainty on effective fi~ nuclear capture has a negligible 
impact on the final measurements. 

Sensitivity to the \i~ content of the data is obtained by 
simultaneously analyzing two samples: those with only a 
muon candidate event, and events consistent with a muon 
and its decay electron. Therefore, this analysis takes as 
signal all and charged-current events. Apart from 
the requirement of cither one or two subevents, the event 
selection for this analysis closely follows that described in 



Section IV B with a few changes appropriate to different 
backgrounds and a higher sensitivity to Michel detection 
efficiency. Table IV details the and charged-current 
purity of the two samples after each cut. 

The primary samples of this analysis are separated 
by Cut 1, where CC events have an enhanced con- 
tribution in the single subevent sample due to pT cap- 
ture. Cuts 2-5 are common to the analysis presented in 
the main body of this work and are motivated in Sec- 
tion IV B. Cuts 6 and 8 reduce the neutral-current back- 
ground in the single-subevent sample: Figure 9 shows 
neutral-current single tt events are largely rejected by 
the requirement on the fi/e log- likelihood variable, while 
Cut 7 further reduces their contribution. Cut 7 uses 
the observed muon kinematics and the stopping power 
of mineral oil for minimum-ionizing particles to calcu- 
late where the muon will stop. This cut removes Michcls 
produced near the optical barrier where Michel detec- 
tion efficiency decreases rapidly with radius and is thus 
sensitive to modeling, while Michel detection is constant 
below 500 cm in this variable. Cut 7 also enhances 
purity due to kinematic differences between and 
CCQE, where the more forward-going nature of the p + 
from P M interactions preferentially stop at high radius 
in the downstream region of the detector. A summary 
of nucleon-level interactions contributing to the selected 
subevent samples is given in Table V. 



Data 

Total MC 

~| Charged current 
| Neutral-current 171 
Other neutral current 




0.15 



FIG. 9: (Color online) The log-likelihood fi/e particle-ID vari- 
able in the single-subevent sample. Events with a muon-like 
score of 0.02 and higher are selected. Expectations are nor- 
malized to flux, and errors shown on data are statistical only. 



3. Calibrations using neutrino-mode data 

Charged-current and events without final-state 
pions typically have two subevents: one from the primary 
(i and another from its decay electron. Two effects deter- 
mine the majority of the migration rate of these events 
from the two-subevent to the one-subevent sample: 
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TABLE IV: Anti-neutrino mode purity in % for all i/ M and charged-current events in the one- and two-subevent samples. A 
pre-cut of generated radius < 550 cm is applied. 



Put ztt 


1 — IT*^tl/"lTl 

UcaLl ip null 


One subevent 


Two subevents 


CC 


v„ CC 


iV CC 


v„ CC 


1 


Subevent cut 


18 


33 


26 


57 


2 


Veto hits < 6 for all subevents 


9 


11 


30 


65 


3 


First subevent in beam window: 4000 < T(ns) < 7000 


9 


11 


29 


65 


4 


Reconstructed vertex radius < 500 cm for first subevent 


8 


11 


29 


65 


5 


Kinetic energy > 200 MeV for first subevent under fi hypothesis 


20 


27 


29 


68 


6 


)x/e log-likelihood ratio > 0.02 for first subevent 


36 


54 


27 


72 


7 


Predicted /i stopping radius < 500 cm 


39 


46 


28 


71 


8 


Q 2 QE > 0.2 GeV 2 


57 


36 


43 


56 



TABLE V: Summary of predicted nucleon-level interactions 
in the anti-neutrino mode subevent samples. The small con- 
tribution from neutral current processes are presented as the 
sum of the v u and v u interactions. 



TABLE VI: A brief description of the neutrino mode subevent 
samples for the same selection described in the previous sec- 



Process 



Contribution (%) to 
One subevent Two subevents 



v^p — > fi + n 
v^n — » fi~p 

i> M (P M )iV -> v^uJN 

other 



31 


49 


48 


36 


3 


5 


7 


7 


1 





3 





4 





3 


3 



nuclear capture and detection efficiency for Michel elec- 
trons. Since an appreciable number (~7%) of charged- 
current events enter the single-subevent sample due to 
Michel detection inefficiencies, the measurement of the 
Un content of the anti-neutrino mode data is sensitive to 
the accuracy of both the Michel detection efficiency and 
the effective fi~ capture rate. The rate of non-detection 
is mostly due to Michel production too close in time with 
the parent muon to be separated by the subevent defini- 
tion. This effect can be seen at low values of the timing 
difference distribution between the two subevents shown 
in Figure 5, Section IV B. 

Fortunately, the neutrino-mode data offers an op- 
portunity to calibrate the migration rate between the 
subevent samples for charged-current events. Due 
to a convolution of flux and cross-section effects [32], 
the neutrino-mode subevent samples are mostly due to 
charged-current interactions. Table VI shows the pre- 
dicted neutrino species and interaction type contribu- 
tions to the neutrino mode subevent samples. With a 
high-purity charged-current sample, the accuracy of 
Michel detection and effective pT capture in simulation 
can be tested. For charged-current events without 
final-state pions [y^ CC), the number of events in the 
neutrino-mode one-subevent (1SE") and two-subevent 
(2SE^) samples are given by: 



tion. 


Process 


Contribution (%) to 
One subevent Two subevents 


All Uy, charged-current 
All 

All neutral current 


95.4 99.0 
0.4 0.7 
4.3 0.3 



TABLE VII: Calibration summary for Michel detection in- 
efficiency (5) and the rate of effective \x~ nuclear capture 
(/3). Note that both processes cannot be simultaneously con- 
strained. 



Process 


data 


MC 


data/MC 


S 


0.073 


0.074 


0.98 


P 


0.071 


0.073 


0.98 



1SE" = i/ M CC x (6 + P(l -S)) (Al) 
2SE" = i^CCx (1- 6-/3(1 -S))+ N£ (A2) 

where N^ (N£) is the neutral-current contribution to the 
1SE (2SE) sample, S is the Michel detection inefficiency 
and {3 is the effective /i - capture rate described previ- 
ously. The rate for Michel non-detection can be solved 
in terms of the effective \i~ capture rate and the small 
neutral-current contribution: 



5 = 



1SE"-Ni 

1SE"+2SE 1 '-(N5;+N£) 



(A3) 



Noting the symmetry in Equations Al and A2 between 
6 and /3, Equation A3 can also express the effective \x~ 
capture rate in terms of Michel detection with S -H- /3. 
Table VII gives values of S and f3 from simulation and 
data based on the observed or predicted event rates in 
the 1SE" and 2SE" samples. 

As the charged-current migration rate to the single- 
subevent sample is due to a convolution of Michel detec- 
tion and effective [i~ capture, the processes cannot be 
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simultaneously calibrated with the neutrino-mode data 
- that is, for example, the calibration of 5 assumes the 
MC valuation of /3 is correct. Future experiments may 
have the ability to separate the two processes by examin- 
ing the low-energy region of the Michel spectrum, where 
the contribution from events following \x~ capture is en- 
hanced. As the calibration results shown in Table VII are 
quite mild and within systematic uncertainties, this pro- 
cedure gives confidence in the ability to unambiguously 
measure the Vn content of the anti-neutrino mode data 
using \jT capture. 

The high-statistics neutrino-mode data also allows for 
a stability check of the ratio of samples one subevent / 
two subevents, and four sequential sample periods are 
consistent within one standard deviation. 



4. Measurement and Systematic Errors 



TABLE VIII: Central-value results for scale factors relative 
to MC expectation for the and charged-current content 
of the anti-neutrino mode data. 



Parameter 


Calibrated 
process 


E$ E range (GeV) 
< 0.9 > 0.9 All 




S 


0.78 


0.79 


0.78 


OLv 


P 


0.78 


0.79 


0.78 




A A T /~\ V* O IT /~\ 


0.78 


0.79 


0.78 




5 


1.16 


1.15 


1.16 




P 


1.16 


1.15 


1.16 




Average 


1.16 


1.15 


1.16 




— Eb) E^ - 


- (E% - 2M p E B + ml - 


h AM 2 ) 




2 [(M p - 


E B ) - E^ + 


cos 9^} 


(A12) 



The flux is measured by adjusting the MC predic- 
tion of the and content to match the data in re- 
gions of reconstructed energy for the subevent samples. 
Following the conventions of Equations Al and A2 and 
introducing CC for the charged-current content, 
the predicted and z/ M contributions to the subevent 
samples in anti-neutrino mode are defined as 



1SE 
'MC 

2SE 
y MC 
-1SE 
y MC 
-2SE 
y MC 



= ^CC x (6 + p(l - 6)) (A4) 

= Uf.CC x (1 - S - 0(1 - 6)) (A5) 

= D^CC x 6 (A6) 

= i/ M CC x (1 - S) (A7) 



Then the single- ("1SE P ") and two-subevent ("2SE S ") 
data samples in anti-neutrino mode are given by 



1SE P = a v x 



Otr, X V 



1SE 
MC 



2SE P 



a v x v 



2SE 
MC 



+ OLy X "MC 



(A8) 
(A9) 



where a v and a p are scale factors for the and 
charged-current content, respectively, to be measured in 
this analysis and the neutral-current content (N£ and Nf) 
include contributions from both and . Equations A8 
and A9 can be solved for a u and a p : 



(1SE P 



iN l7^MC 



(2SE P 



.-.2SE..1SE 
"MC "MC 

(lSE^-ND^c 5 



1SE,.2SE 
MC ^MC 



V_MC (A1Q) 



2SE r .lSE 
MC "MC 



(2SE^-N^)^ (An) 



v- 



1SE.-.2SE 
MC "MC 



where, Eb = 30 MeV is the binding energy, m M is the 
muon mass, AM 2 = Af 2 - A/ 2 , where M n (M p ) is the 
neutron (proton) mass, is the muon momentum, and 
9^ is the outgoing muon angle relative to the incoming 
neutrino beam. This reconstruction assumes CCQE 
interactions with at-rest, independently-acting nucleons. 
Though this is a model-dependent valuation of the neu- 
trino energy, complicated further by the significant non- 
Vfj, CCQE content, separating the samples into exclusive 
regions of E® E nevertheless affords statistical sensitiv- 
ity to the accuracy of the simulated flux spectrum. The 
three energy regions explored arc E® E < 0.9 GeV, 
Ej E > 0.9 GeV, and an inclusive sample. The statis- 
tics of the single-subevent sample prohibit the analysis 
of more than two exclusive E$ E regions. As described in 
the previous section, the calibration from the neutrino- 
mode data is ambiguous between Michel detection and 
the effective [i~ capture model. As these effects change 
the expectations for ^mc^mc^mc anc ^ ^mc m differ- 
ent ways, the measurement of a v and a t is, in principle, 
sensitive to which rate is calibrated. In the absence of a 
compelling reason to choose one over the other, the final 
evaluations for a v and ap are taken to be the average of 
the two calculations assuming each rate is calibrated. A 
calibration uncertainty spanning the difference in the two 
measurements is added to the systematic errors discussed 
next. The central values for a v and otp are presented in 
Table VIII. 

Systematic uncertainties on a v and a e are evaluated 
by assigning relevant errors to the physics processes con- 
tributing to the subevent samples and observing how the 
measurement changes as the channels arc varied within 
their uncertainty These uncertainties are treated as un- 
correlatcd, so the uncertainty on a v , for example, due to 
physics processes P\ , ■ ■ ■ , Pjv is simply 



To check the modeling of the flux spectrum, this 
measurement is performed in three regions of recon- 
structed energy E® E , defined as 



N / a \ 2 



(A13) 
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Tabic IX shows the errors assigned to the various 
contributing processes and their propagated uncertainty 
onto <x v and a p . The most important process for extract- 
ing the flux measurement is the CCQE interaction, 
and its cross section and assigned uncertainty reflect the 
measurement and accuracy of the MiniBooNE result [19]. 
The same is true for the and neutral-current sin- 
gle 7T° channels [72]; however the error is increased to 
recognize a possible rate difference in these interactions 
between the cross-section measurements and this analysis 
due to using the opposite side of the log-likelihood vari- 
able shown in Figure 9. The and charged-current 
single charged 7r channels are adjusted to reflect the 
measurement [19] and their uncertainty is increased to 
recognize the extrapolation to the processes. Treat- 
ing the uncertainties on the processes constrained by 
MiniBooNE data as uncorrelated ignores a common de- 
pendence on the neutrino-mode flux uncertainties and 
a small cancelation of errors that could be propagated 
onto a v and ap is ignored. The neutral-current elas- 
tic process is also constrained by MiniBooNE data [73], 
while the neutral-current charged-pion production pro- 
cesses are completely unconstrained and so the assigned 
uncertainty is large. Preliminary results for the CCQE 
process [36] informs the choice of a 20% uncertainty rela- 
tive to the RFC model with M A = 1.35 GeV. With these 
systematic uncertainty assumptions, as seen in Table IX, 
the uncertainty on the main result of this work a v is 
dominated by statistics and the CCQE cross section. 
As the CCQE process is directly constrained by Mini- 
BooNE data, the measurement of the flux scale a v 
features negligible model dependence. Table X summa- 
rizes the measurements of a v and a e . 

As the cross sections for the dominant processes 
have been applied to simulation, the deviation from unity 
for a v represents the accuracy of the highly-uncertain 
flux prediction in anti-neutrino mode. As the bulk 
of the flux prediction is constrained by the HARP 
data, the a p scale factor is representative of the level of 
cross-section agreement between data and the RFC with 
M A = 1.35 GeV for the ^ CCQE process. 



5. Summary 

This Appendix presents a measurement of the flux 
in anti-neutrino mode using a non-magnetized detector. 
The results arc consistent with and complementary to 
the two measurements in Ref. [32]. A summary of the 
results from all three analyses is shown in Figure 3, Sec- 
tion IV A. As no energy dependence among the mea- 
surements is observed, the simulation of the flux in 
anti-neutrino mode, which is unconstrained by the HARP 
hadroproduction data, appears to be roughly 20% high in 
normalization, while the flux spectrum is well-modeled. 

These techniques could also aid future neutrino exper- 
iments that will test for CP violation in the lepton sector 
using large unmagnctized detectors such as NO^A [5], 



T2K [4], LBNE [7], LACUNA [8], and Hyper-K [9]. In 
particular, the precision of ~ 15% in the determination 
of the flux of the anti-neutrino mode beam using /i~ 
capture obtained here could easily be surpassed and the 
flux spectrum more rigorously checked by future exper- 
iments housing heavier nuclei. As an example, the fi~ 
capture rate on 40 Ar exceeds 70% [65], almost afford- 
ing event-by-event discrimination of the fi charge without 
a magnetic field. Detector-specific complications arising 
from 7r//x identification and Michel detection should not 
reduce sensitivity to the ji charge dramatically. 



Appendix B: Model-dependent measurements for 
interactions on CH2 

This Appendix presents MiniBooNE CCQE cross- 
section measurements that are explicitly dependent on 
CCQE interaction assumptions. These measurements 
include all CCQE interactions as signal, while Ap- 
pendix C gives cross sections treating the hydrogen 
CCQE component as background. All results are tab- 
ulated in Appendix D. 



1. Total cross section 

As the energy distribution of the incident beam is 
quite broad (Figure 1), the a priori knowledge of the 
neutrino energy is highly uncertain on an event- by- event 
basis. If hadronic reconstruction is unavailable, it is typ- 
ical for neutrino experiments to reconstruct the neutrino 
energy of events in the CCQE sample assuming scat- 
tering off of at-rest and independently-acting nucleons 
("EQ e ") based solely on the outgoing lepton kinemat- 
ics (Equation A12). Finding the neutrino energy in this 
way is often used to measure neutrino oscillation param- 
eters, in particular the mass splitting, and it has been 
argued elsewhere that the assumptions implicit in this 
reconstruction significantly bias these measurements due 
to ignored nuclear effects [10, 11]. 

Apart from the bias in the reconstructed energy dis- 
tribution, a measurement of the absolute cross section 
over the observed energy range additionally suffers from 
model dependence through the unfolding procedure. The 
total cross section is typically computed by unfolding the 
reconstructed neutrino energy to the "true" energy dis- 
tribution, and this correction is dependent on both the 
nuclear model used and detector resolution effects. This 
is the main reason MiniBooNE has generally opted to re- 
port cross sections in terms of observed kinematics. Due 
to these measurement biases, the MiniBooNE CCQE 
absolute cross section is not the main result of the work 
but is provided here for historical comparisons. 

A consequence of the unfolding bias is that one should 
exercise caution in comparing theoretical calculations to 
these results. A strict comparison with these data and an 
external model involves finding the total cross section as 
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TABLE IX: Uncertainty summary for this analysis. Included are the assumed error on physics processes and their contribution 
to the total error in a v and in the regions of reconstructed neutrino energy studied. The statistics of the v-mode data enter 
the uncertainty from the calibration procedure described in Section IV B. 



uU L1I Ut; 


Fractional 
uncertainty (%) 


Uncertainty contribution to a 
E$ E < 0.9 GeV E$ E > 0.9 GeV 


V 

All 


Uncertainty contribution to a 
E$ E < 0.9 GeV E$ E > 0.9 GeV 


All 


u^n — y n~p 


10 


0.07 


0.08 


0.07 


0.00 


0.00 


0.00 


Vy,p — > fi + n 


20 


0.04 


0.02 


0.03 


0.20 


0.20 


0.21 




) 20 


0.04 


0.05 


0.04 


0.02 


0.02 


0.01 




30 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


v^)N —¥ v»{p,j)Nir° 


25 


0.02 


0.01 


0.01 


0.01 


0.01 


0.01 




50 


0.05 


0.02 


0.01 


0.03 


0.03 


0.01 


H~ capture 


2.8 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


i^-mode statistics 




0.10 


0.11 


0.08 


0.08 


0.08 


0.06 


z^-mode statistics 




0.04 


0.05 


0.04 


0.03 


0.03 


0.03 


All 


0.14 


0.16 


0.12 


0.22 


0.22 


0.22 



TABLE X: Summary of measurements for the flux scale 



a v and the 


rate scale ap. 






Parameter 


< 0.9 


E$ E range (GeV) 
> 0.9 


All 


u v 


0.78 ± 0.14 


0.79 ± 0.16 


0.78 ± 0.12 




1.16 ± 0.22 


1.15 ± 0.22 


1.16 ± 0.22 



a function of E® E using the generated muon kinematics, 
and subsequently unfolding this distribution according to 
the RFC An example of this procedure can be found in 
Ref. [74]. 

The flux-unfolded CCQE cross section per nucleon 



is calculated assuming: 



x10" 



0.4 



MiniBooNE v„ CCQE data (CHJ 

Free nucleon, M' 4 " = 1.02 GeV 

Mf c (M*™) = 1.35 (1.02) GeV,K = 1.007 

M? C = M*-: 



= 1.02 GeV,K= 1.000 




0.6 



0.8 



1 1.2 

E OE,RFG (QeV) 



1.4 



1.6 
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where the same conventions used in Equation 3 apply 
here with a few exceptions: as mentioned, the unfolding 
matrix Uij here connects the reconstructed neutrino en- 
ergy (inferred from the observed /i kinematics via Equa- 
tion A12) to the generated distribution, and the flux term 
<i>i refers to the flux exclusive to the i th neutrino en- 
ergy bin. Figure 10 compares the observed total cross 
section to a few predictions from the RFC 



2. Momentum transfer 

Another important quantity for CCQE interactions is 
the squared four- momentum transfer Q 2 = (p v — p^) . 
However, again ignorance of the incoming neutrino en- 
ergy prevents a clean measurement of this variable. As 
in the case for E® E (Equation A12), if only lepton kine- 
matics are available the distribution can be inferred by 
assuming CCQE scattering with an at-rest, independent 
nucleon: 



FIG. 10: (Color online) Per-nucleon total cross section for 
MiniBooNE CCQE data including the hydrogen scatter- 
ing component. The distribution is labeled Eg ' to rec- 
ognize the dependence on the assumptions inherent in both 
the reconstruction and in the unfolding model. Total errors 
are shown with data. Numerical values are provided in Ta- 
ble XVI. 



l2 2 



2E$ E (E^-p fl cose fi ) (B2) 



where E^ , p M and m M refer to the muon energy, momen- 
tum and mass, respectively. The value of the axial mass is 
typically extracted from the shape of this distribution, so 
the differential cross section with respect to this variable 
is provided for historical comparisons despite the recon- 
struction assumptions. However, to minimize the model 
dependence of this cross-section configuration, the recon- 
structed distribution of Qq E is corrected to true Qqe ~ 
that is, Equation B2 with the generated muon kinemat- 
ics. In this way, the unfolding procedure only corrects for 
muon resolution effects and is not biased by the CCQE 
interaction model. Note that truth-level Qqe is omv 
the same as the squared four-momentum transfer up to 
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the naive reconstruction assumptions. This choice is not 
typically made and so comparisons with similar cross sec- 
tions from other experiments should be made with care. 

The flux-folded, single-differential cross section 
da /dQq E calculated in the same manner as the double- 
differential cross section (Equation 3) but for a single 
dimension: 



da 



v. Uij (dj - bj) 



dQ 2 



QE 



(AQ 2 



QE 



)iei*JV" 



(B3) 



where the same conventions used in Equation 3 apply. 
Figure 11 compares the results with shape uncertainty 
to predictions from the RFG normalized to data. 
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FIG. 11: (Color online) Per-nucleon single-differential cross 
section for MiniBooNE CCQE data including hydrogen 
CCQE events. The RFG predictions are normalized to the 
observed total cross section J j^—dQg E and the relative 

scales are indicated. All predictions assume an effective axial 
mass of 1.026 GeV for the hydrogen scattering component. 
Shape errors are shown with data. Numerical values are pro- 
vided in Table XV. 

The conventions used to calculate <t(E„) and -tM— are 

the same used to calculate the corresponding CCQE 
cross sections reported in Rcf. [19]. 

Appendix C: Model-dependent measurements for 
CCQE interactions on carbon 

Following the same definitions for the total and single- 
differential (Equations Bl, and B3, respectively) cross 
sections, results for CCQE on carbon are obtained af- 
ter an RFG prediction for the hydrogen component 
with Ma = 1.026 GeV. In this configuration, the total 
per-nucleon cross section for CCQE interactions on car- 
bon from both the MiniBooNE and analyses may 
be compared to the corresponding NOMAD results, and 
this is shown in Figure 12. 



FIG. 12: (Color online) Total i/ M and CCQE cross sec- 
tions for the MiniBooNE and NOMAD experiments, shown 
with two choices for the axial mass in the RFG for compar- 
ison. The hydrogen content has been subtracted from the 
MiniBooNE data, and total uncertainties are shown. Note 
the two experiments use difference detector technologies and 
so naturally assume different topologies in defining CCQE 
events. Therefore, conclusions should be drawn with care. 
NOMAD data taken from Ref. [23], and MiniBooNE data 
taken from Ref. [19]. Numerical values for the MiniBooNE 
cross section are provided in Table XXI. 



Appendix D: Tabulation of results 

This Appendix provides numerical values for the ob- 
served flux and all cross-section results presented in 
this work. In addition, each cross section is accompa- 
nied by a measurement of the "CCQE-likc" background. 
This corresponds to the amount of processes that con- 
tribute to the analysis sample, mostly CCl7r~, that have 
been subtracted in the main analysis. Note that in or- 
der to facilitate comparisons with the predictions of anti- 
neutrino CCQE and CCQE-likc processes, the CCQE- 
like measurements exclude the content of the sub- 
tracted data. The CCQE-likc background is calculated 
for the various cross sections (Equations 3, Bl, and B2) 
by replacing (dj — bj) with the subset of bj correspond- 
ing to only those induced by v^. Note also these mea- 
surements are normalized to the total number of proton 
targets in the detector, even though the dominant inter- 
action of CC17T - has nucleon-level interactions with neu- 
trons as well. This configuration is chosen for consistency 
with the Vy, CCQE-likc background measurements, which 
were normalized to the number of neutron targets in the 
Vy CCQE analysis [19]. As the CCQE-like cross sections 
on mineral oil and carbon differ only by the inclusion of 
the hydrogen content, the amount of Vp hydrogen CCQE 
subtracted from the data (in the case of the latter calcu- 
lation) can be found by taking the difference of these two 
cross sections. To find the calculated per-nucleon hy- 
drogen CCQE cross section, this difference should also be 
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scaled by the ratio of total protons targets to quasi-free 
proton targets, 2.03 x 10 32 /0.70 x 10 32 = 2.9. 

1. Anti-neutrino flux 

Section II A describes the flux prediction, and Ta- 
ble XI lists the predicted flux per POT in 50 MeV 
wide bins of energy up to 3 GcV. These values normal- 
ized to the observed exposure of 10.1 x 10 20 POT are 
shown in Figure 1. 

2. Cross-section results on CH2 

All measurements in this section include the quasi-free 
hydrogen CCQE scattering component and is therefore 
less model-dependent compared to the results given in 
Appendix D 3, where the RFG model is relied on to sub- 
tract their contribution. Shape uncertainties are pro- 
vided for the double- and single-differential cross-section 
measurements, and these values should be used along 
with the total normalization uncertainty of 12.9% in the 
context of a fit to these distributions. The total uncer- 
tainty, including errors affecting both shape and normal- 
ization, is provided for the total cross section. 

Numerical values for the MiniBooNE P M CCQE cross 
section including the hydrogen content is given in Ta- 
ble XII, while Table XIII provides the uncertainty on 



the shape of these data. These tables correspond to Fig- 
ures 7 and 8. The CCQE-like background is reported in 
Table XIV. 

The single-differential cross section -tM— measurement 

6 d Q QE 

with shape uncertainty and CCQE-like background is 
given in Table XV. 

The total cross section including the hydrogen content 
is given in Table XVI. As discussed in Appendix B, these 
results are dependent on both the E^ E reconstruction 
assumptions (Equation A12) and the nuclear model used. 
To recognize these dependencies, the neutrino energy is 
labeled here as "£QE,Rfg» _ 



3. Cross-section results on C 

Tables in this section report the cross-section results 
reliant on the RFG subtract hydrogen CCQE events from 
the data. In addition to the shape errors provided with 
each measurement, a normalization uncertainty of 17.2% 
is applicable here. 

The MiniBooNE CCQE double-differential cross 
section, shape uncertainty, and CCQE-like cross section 
treating the hydrogen CCQE content as background is 
given in Tables XVII, XVIII, and XIX, respectively. Ta- 
bles XX and XXI provide the same information for the 
single-differential and total cross section. 



21 



E v bin flux 
(GeV) (^/POT/GeV/cm 2 ; 



E v bin flux 

(GeV) (pn /POT/50 MeV/cm 2 ) 



0.00-0.05 


2.157 


X 


1(T 


12 


1 


00- 


1 


05 


7.658 


X 


io- 


12 


2 


00- 


2 


05 


2.577 


X 


10" 


13 


0.05-0.10 


7.840 


X 


10" 


12 


1 


05- 


1 


10 


6.907 


X 


io- 


12 


2 


05- 


2 


10 


2.066 


X 


10- 


13 


0.10-0.15 


9.731 


X 


10" 


12 


1 


10- 


1 


15 


6.180 


X 


io- 


12 


2 


10- 


2 


15 


1.665 


X 


io- 


13 


0.15-0.20 


1.141 


X 


10" 


11 


1 


15- 


1 


20 


5.505 


X 


io- 


12 


2 


15- 


2 


20 


1.346 


X 


10- 


13 


0.20-0.25 


1.319 


X 


10" 


11 


1 


20- 


1 


25 


4.877 


X 


io- 


12 


2 


20- 


2 


25 


1.081 


X 


10~ 


13 


0.25-0.30 


1.438 


X 


10" 


11 


1 


25- 


1 


30 


4.269 


X 


io- 


12 


2 


25- 


2 


30 


8.837 


X 


io- 


14 


0.30-0.35 


1.477 


X 


10" 


11 


1 


30- 


1 


35 


3.686 


X 


io- 


12 


2 


30- 


2 


35 


7.136 


X 


io- 


14 


0.35-0.40 


1.479 


X 


10" 


11 


1 


35- 


1 


40 


3.151 


X 


io- 


12 


2 


35- 


2 


40 


5.707 


X 


10- 


14 


0.40-0.45 


1.500 


X 


10" 


11 


1 


40- 


1 


45 


2.678 


X 


io- 


12 


2 


40- 


2 


45 


4.620 


X 


10- 


14 


0.45-0.50 


1.485 


X 


10" 


11 


1 


45- 


1 


50 


2.262 


X 


io- 


12 


2 


45- 


2 


50 


3.778 


X 


io- 


14 


0.50-0.55 


1.447 


X 


10" 


11 


1 


50- 


1 


55 


1.898 


X 


io- 


12 


2 


50- 


2 


55 


3.028 


X 


10" 


14 


0.55-0.60 


1.406 


X 


10" 


11 


1 


55- 


1 


60 


1.580 


X 


io- 


12 


2 


55- 


2 


60 


2.412 


X 


io- 


14 


0.60-0.65 


1.345 


X 


10" 


11 


1 


60- 


1 


65 


1.311 


X 


io- 


12 


2 


60- 


2 


65 


1.977 


X 


10" 


14 


0.65-0.70 


1.287 


X 


10" 


11 


1 


65- 


1 


70 


1.083 


X 


io- 


12 


2 


65- 


2 


70 


1.638 


X 


io- 


14 


0.70-0.75 


1.221 


X 


10" 


11 


1 


70- 


1 


75 


8.917 


X 


io- 


13 


2 


70- 


2 


75 


1.323 


X 


10- 


14 


0.75-0.80 


1.152 


X 


10" 


11 


1 


75- 


1 


80 


7.285 


X 


io- 


13 


2 


75- 


2 


80 


1.038 


X 


10" 


14 


0.80-0.85 


1.075 


X 


10" 


11 


1 


80- 


1 


85 


5.941 


X 


io- 


13 


2 


80- 


2 


85 


8.707 


X 


10" 


15 


0.85-0.90 


9.980 


X 


io- 


12 


1 


85- 


1 


90 


4.834 


X 


io- 


13 


2 


85- 


2 


90 


6.981 


X 


io- 


15 


0.90-0.95 


9.177 


X 


10" 


12 


1 


90- 


1 


95 


3.937 


X 


io- 


13 


2 


90- 


2 


95 


6.078 


X 


10" 


15 


0.95-1.00 


8.411 


X 


10" 


12 


1 


95- 


2 


00 


3.180 


X 


10" 


13 


2 


95- 


3 


00 


5.111 


X 


io- 


15 



E v bin Vfj, flux 

(GeV) (^/POT/GeV/cm 2 ) 



TABLE XI: Predicted Vu flux at the MiniBooNE detector in anti-neutrino mode. 



TABLE XII: The MiniBooNE z/ M CCQE double-differential cross section on mineral oil in units of 10 41 cm 2 /GeV. Data is 
given in 0.1 GeV bins of T M (columns) and 0.1 bins of cos# M (rows). 
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TABLE XIII: Shape uncertainty in units of 10 41 cm 2 /GeV on the MiniBooNE CCQE double-differential cross section on 
mineral oil. The total normalization uncertainty is 12.9%. 
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TABLE XIV: CCQE-like background in units of 10 41 cm 2 /GeV to the MiniBooNE CCQE double-differential cross section 
on mineral oil. In this configuration, the hydrogen scattering component is treated as signal and is not included in the CCQE-like 
background. 
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TABLE XV: The MiniBooNE v„ CCQE single differential cross section d ° on mineral oil, shape error, and CCQE-like 
background in units of cm 2 /GeV 2 . The total normalization error is 12.9%. 
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TABLE XVI: The MiniBooNE CCQE total cross section on mineral oil, errors, and predicted CCQE-like background in 
bins of £? E ' RFG 



ggETHTO (QeV) a (cm 2 ) shape error (cm 2 ) total error (cm 2 ) CCQE-like background (cm 2 )" 
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10" 


-39 


3 


132 


X 


10" 


-40 


3.132 


X 


10" 


-40 


7.205 


X 


10" 


-40 


0.60-0.65 


2.542 


X 


10" 


-39 


2, 


,888 


X 


10" 


-40 


3.086 


X 


10" 


-40 


7.534 


X 


10" 


-40 


0.65-0.70 


2.753 


X 


10" 


-39 


3 


,279 


X 


10" 


-40 


3.363 


X 


10" 


-40 


7.799 


X 


10" 


-40 


0.70-0.75 


2.932 


X 


10" 


-39 


2, 


736 


X 


10" 


-40 


3.213 


X 


10" 


-40 


8.048 


X 


10" 


-40 


0.75-0.80 


3.098 


X 


10" 


-39 


2, 


,895 


X 


10" 


-40 


3.399 


X 


10" 


-40 


8.334 


X 


10" 


-40 


0.80-0.90 


3.374 


X 


10" 


-39 


2, 


691 


X 


10" 


-40 


3.596 


X 


10" 


-40 


8.848 


X 


10" 


-40 


0.90-1.00 


3.780 


X 


10" 


-39 


2, 


,385 


X 


10" 


-40 


4.165 


X 


10" 


-40 


9.517 


X 


10" 


-40 


1.00-1.10 


4.171 


X 


10" 


-39 


2, 


433 


X 


10" 


-40 


4.970 


X 


10" 


-40 


1.017 


X 


10" 


-39 


1.10-1.30 


4.631 


X 


10" 


-39 


3 


,755 


X 


10" 


-40 


7.339 


X 


10" 


-40 


1.087 


X 


10" 


-39 


1.30-1.50 


5.510 


X 


10" 


-39 


8, 


,576 


X 


10" 


-40 


1.297 


X 


10" 


-40 


1.261 


X 


10" 


-39 


1.50-2.00 


6.654 


X 


10" 


-39 


1 


910 


X 


10" 


-39 


2.402 


X 


10" 


-40 


1.388 


X 


10" 


-39 
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TABLE XVII: The MiniBooNE z? M CCQE double-differential cross section on carbon. The units are 1(T 41 cm 2 /GeV. 



a T.AGeV) 


0.2,0.3 0.3,0.4 0.4,0.5 0.5,0.6 0.6,0.7 0.7,0.8 0.8,0.9 0.9,1.0 1.0,1.1 1.1,1.2 1.2,1.3 1.3,1.4 1.4,1.5 1.5,1.6 1.6,1.7 1.7,1.8 1.8,1.9 1.9,2.0 


+0.9, + 1.0 


213.4 


311.5 


517.6 


710.8 


848.8 


969.0 1015 998.8 973.0 875.0 726.3 638.5 547.9 484.0 


+0.8, +0.9 


301.1 


448.9 


664.0 


800.9 


839.5 


743.8 667.6 518.4 389.7 238.6 126.7 92.13 


+0.7, +0.8 


304.7 


406.6 


535.9 


521.4 


456.9 


355.0 233.2 135.1 61.71 


+0.6, +0.7 


295.3 


343.1 


362.1 


322.4 


228.8 


129.0 59.01 — 


+0.5, +0.6 


294.6 


299.5 


238.3 


178.4 


109.0 


47.10 — — — — — — — — — — — — 


+0.4, +0.5 


228.6 


214.0 


145.7 


88.79 


34.79 




+0.3, +0.4 


178.9 


143.2 


94.80 


33.00 






+0.2, +0.3 


156.9 


83.82 


34.40 


11.86 






+0.1, +0.2 


103.6 


52.63 


17.23 








0.0, +0.1 


89.47 


47.39 


8.768 








-0.1,0.0 


80.99 


19.33 










-0.2,-0.1 


53.95 


12.83 










-0.3,-0.2 


30.79 












-0.4,-0.3 


24.58 












-0.5,-0.4 


27.97 












-0.6,-0.5 


24.58 












-0.7,-0.6 














-0.8,-0.7 














-0.9,-0.8 














-1.0,-0.9 















TABLE XVIII: Shape uncertainty for the CCQE double-differential cross section on carbon. The units are 10 41 cm 2 /GeV, 
and the total normalization uncertainty is 17.2%. 





0.2,0.3 0.3,0.4 0.4,0.5 0.5,0.6 0.6,0.7 0.7,0.8 0.8,0.9 0.9,1.0 1.0,1.1 1.1,1.2 1.2,1.3 1.3,1.4 1.4,1.5 1.5,1.6 1.6,1.7 1.7,1.8 


i 1.8,1.9 1.9,2.0 


+0.9, + 1.0 


97.99 


100.4 


123.6 


137.2 


132.8 


124.2 118.4 116.8 133.3 145.9 152.7 173.4 193.7 308.1 




+0.8, +0.9 


89.94 


97.56 


105.6 


103.8 


111.1 


78.41 70.51 69.88 75.89 82.72 57.19 58.51 — — 




+0.7, +0.8 


77.84 


73.26 


77.15 


68.58 


58.20 


49.52 46.88 40.61 34.87 — — — 




+0.6, +0.7 


71.35 


57.41 


53.55 


45.39 


39.75 


31.00 36.81 — 




+0.5, +0.6 


65.99 


48.22 


37.81 


32.38 


28.76 


22.26 — — — — — — — — — — 




+0.4, +0.5 


50.11 


38.05 


29.92 


30.11 


17.85 






+0.3, +0.4 


42.73 


30.20 


24.09 


28.22 








+0.2, +0.3 


39.17 


22.70 


19.03 


11.05 








+0.1, +0.2 


32.42 


17.09 


13.77 










0.0, +0.1 


28.41 


18.48 


10.54 










-0.1,0.0 


26.47 


12.27 












-0.2,-0.1 


21.54 


10.22 












-0.3,-0.2 


16.94 














-0.4,-0.3 


22.87 














-0.5,-0.4 


19.70 














-0.6,-0.5 


24.78 














-0.7,-0.6 
















-0.8,-0.7 
















-0.9,-0.8 
















-1.0,-0.9 
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TABLE XIX: The QE-like background subtracted from the double-differential cross section on carbon. The z/ M CCQE 
interactions with hydrogen are treated as background in this calculation, and so their contribution is included here. The units 
are 1(T 41 cm 2 /GeV. 





0.2,0.3 0.3,0.4 0.4,0.5 0.5,0.6 0.6,0.7 0.7,0.8 0.8,0.9 0.9,1.0 1.0,1.1 1.1,1.2 1.2,1.3 1.3,1.4 1.4,1.5 1.5,1.6 1.6,1.7 1.7,1.8 


i 1.8,1.9 1.9,2.0 


+0.9, + 1.0 


222 1 


1 no s 


O O O . )) 


DDU.O 


7^R O 7fifi fl 7^9 7 fiOQ Q fif| 1 R /IQfi K ^Q9 9 ^n7 1 997 9 1 ^7 1 

(oo.y (DD.u ioz.i ooy.o dui.d lyo.o oy^.z ou/ .1 zz/.z 10/. 1 — — 




+0.8, +0.9 


1 R7 7 


Tl /I Q 




^sin i 

OOU. 1 


■3^7(1 9QQD 9/1 n A 189 118/1 7/1 /II /in 

oo i .u zyy .u z^u.^t ioz.o 110.4 14.00 *± 1 . iu — — — — — 




+0.7, +0.8 


142 5 


211 1 


217 


190 8 


1 ^1 s 1 s fi^ 1 ^ ^7 sn 

IOI.O lUO.o DO.IO O/.oU 




+0.6, +0.7 


107 9 


144 5 


128 6 


94 78 


fi9 ^ Qn 1 7 SS 
oz.zo oo.yu 1 1 .00 




+0.5, +0.6 


oz.Oo 






46 40 


no 

zo. yo 




+0.4, +0.5 


R9 Q9 


R9 SS 
Oz.Oo 


41 1 1 


99 9 9 






+0.3, +0.4 


47.91 


41.44 


23.83 








+0.2, +0.3 


33.44 


25.40 


11.47 








+0.1, +0.2 


23.84 


16.09 










0.0, +0.1 


18.46 


10.31 










-0.1,0.0 


13.22 












-0.2,-0.1 














-0.3,-0.2 














-0.4,-0.3 














-0.5,-0.4 














-0.6,-0.5 














-0.7,-0.6 














-0.8,-0.7 














-0.9,-0.8 














-1.0,-0.9 















TABLE XX: The MiniBooNE Vn CCQE single differential cross section ' %f on carbon, shape error, and CCQE-like back- 
ground in units of cm 2 /GeV 2 . The CCQE content is treated as background, and the total normalization error is 17.2%. 



Qqe (GeV 2 ) dcp- E shape uncertainty CCQE-like background 



0.00- 


0.05 


6.076 


X 


10" 


39 


1.002 


X 


10" 


-40 


7.809 


X 


10" 


39 


0.05- 


0.10 


8.769 


X 


10" 


39 


5.286 


X 


nr 


-40 


5.774 


X 


10" 


:i<) 


0.10- 


0.15 


7.495 


X 


10" 


-39 


2.846 


X 


10" 


-41) 


3.727 


X 


10" 


-39 


0.15- 


0.20 


6.202 


X 


10" 


39 


2.384 


X 


Kr 


-40 


2.569 


X 


10" 


39 


0.20- 


0.25 


4.921 


X 


10" 


39 


2.160 


X 


10" 


-41) 


1.812 


X 


10" 


-39 


0.25- 


0.30 


3.691 


X 


10" 


39 


1.915 


X 


nr 


-40 


1.284 


X 


10" 


-39 


0.30- 


0.35 


2.838 


X 


10" 


39 


2.079 


X 


10" 


-40 


9.322 


X 


10" 


40 


0.35- 


0.40 


2.297 


X 


10" 


-39 


2.265 


X 


10" 


-40 


6.851 


X 


10" 


40 


0.40- 


0.45 


1.786 


X 


10" 


-39 


1.676 


X 


10" 


-40 


5.077 


X 


10" 


40 


0.45- 


0.50 


1.418 


X 


10" 


39 


1.912 


X 


10" 


-40 


3.841 


X 


10" 


40 


0.50- 


0.60 


9.995 


X 


10" 


4(1 


1.536 


X 


10" 


-40 


2.603 


X 


10" 


40 


0.60- 


0.70 


5.981 


X 


10" 


-40 


1.142 


X 


lo- 


-40 


1.557 


X 


10" 


40 


0.70- 


0.80 


3.757 


X 


10" 


4(1 


1.217 


X 


ur 


-40 


9.733 


X 


10" 


-41 


0.80- 


1.00 


2.104 


X 


10" 


-40 


8.737 


X 


nr 


-41 


5.225 


X 


10" 


41 


1.00- 


1.20 


9.186 


X 


10" 


41 


6.992 


X 


10" 


-41 


2.382 


X 


10" 


41 


1.20- 


1.50 


3.099 


X 


10" 


41 


4.172 


X 


10- 


-41 


1.014 


X 


10" 


41 


1.50- 


2.00 


7.414 


X 


10" 


42 


2.004 


X 


10" 


-41 


3.167 


X 


10" 


4_> 
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TABLE XXI: The MiniBooNE CCQE total cross section on carbon, errors, and predicted CCQE-like background in bins of 

£;QE,RFG 



E ya,x,r<j (G e v) cr (cm 2 ) shape error (cm 2 ) total error (cm 2 ) CCQE-like background (cm 2 ) 






.40-0. 


,45 


1. 


.808 


X 


10" 


-M 


6. 


211 


X 


10" 


-40 


6, 


211 


X 


10" 


-4U 


1 


.127 


X 


10" 







.45-0 


.50 


1 


.890 


X 


10" 


-39 


4. 


,372 


X 


ur 


-40 


4 


372 


X 


10" 


-40 


1 


.224 


X 


10" 


-:«> 





.50-0 


.55 


2. 


019 


X 


10" 


-39 


4. 


319 


X 


10" 


-40 


4, 


319 


X 


10- 


-40 


1 


.309 


X 


io- 


-39 





.55-0 


.60 


2. 


,258 


X 


io- 


-39 


4. 


,068 


X 


10- 


-40 


4 


267 


X 


10- 


-40 


1 


.386 


X 


io- 


-39 





.60-0 


.65 


2. 


501 


X 


10" 


-39 


3. 


,737 


X 


10" 


-41) 


4, 


228 


X 


iir 


-40 


1 


.454 


X 


io- 


-39 





.65-0 


70 


2. 


,728 


X 


io- 


-39 


4. 


,202 


X 


nr 


-40 


4 


483 


X 


10- 


-40 


1 


.512 


X 


io- 


-: J ,9 





.70-0 


.75 


2. 


932 


X 


10" 


-39 


3. 


,528 


X 


10" 


-40 


4, 


335 


X 


iir 


-40 


1 


.575 


X 


io- 


-39 





.75-0 


.80 


3. 


091 


X 


io- 


-39 


3. 


,573 


X 


10- 


-40 


4 


524 


X 


10- 


-40 


1 


.645 


X 


io- 


-39 





.80-0 


.90 


3. 


372 


X 


10" 


-39 


3. 


,368 


X 


10" 


-40 


4, 


806 


X 


10- 


-40 


1 


.753 


X 


io- 


-39 





.90-1 


.00 


3. 


,815 


X 


io- 


-39 


3. 


141 


X 


10- 


-40 


5, 


659 


X 


10" 


-40 


1 


.895 


X 


10" 


-: J ,9 


1 


.00-1 


.10 


4. 


,254 


X 


10" 


-39 


3. 


,250 


X 


10" 


-40 


6, 


702 


X 


lir 


-40 


2 


.022 


X 


10- 


-39 


1 


.10-1 


.30 


4. 


,789 


X 


io- 


-39 


5. 


,144 


X 


10- 


-40 


9. 


.829 


X 


10- 


-40 


2 


.121 


X 


io- 


-39 


1 


.30-1 


.50 


5. 


,784 


X 


10" 


-39 


1. 


159 


X 


io- 


-39 


1. 


741 


X 


lir 


-39 


2 


.378 


X 


10" 


-39 


1 


.50-2 


.00 


7 


086 


X 


io- 


-: J ,9 


2, 


,439 


X 


10- 


-39 


3. 


119 


X 


10 - 


-39 


2 


.482 


X 


io- 


-: J ,9 



